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^ , Central colhsions induced by high energy radioactive beams can be used as a novel 

' means to obtain crucial information about the high density (HD) behaviour of nuclear 
■ symmetry energy. This information is critical for understanding several key issues in astro- 
^ physics. Within an isospin-dependent hadronic transport model using phenomenological 
O ! equations of state (EOS) for dense neutron-rich matter, we investigate several experimen- 
2 I tal probes of the HD behavior of nuclear symmetry energy, such as, the tt~ to tt"'" ratio, 
neutron-proton differential flow and its excitation function. Measurements of these ob- 
servables will provide the first terrestrial data to constrain stringently the HD behaviour 
of nuclear symmetry energy and thus also the EOS of dense neutron-rich matter. 

1. INTRODUCTION 



^ ! Nuclear reactions using rare isotopes has opened up several new frontiers in nuclear 
^ [ sciences [1-3]. In particular, high energy heavy rare isotopes to be available at the future 
Rare Isotope Accelerator (RIA) and the new GSI accelerator facility provide a unique 
opportunity to explore novel properties of dense neutron-rich matter that was not in 
reach in terrestrial laboratories before. This exploration will reveal crucial information 
about the EOS of neutron-rich matter. To understand the latter and its astrophysical 
implications, such as, the origin of elements, structure of rare isotopes and properties of 
neutron stars are presently among the most important goals of nuclear sciences. The EOS 
of neutron-rich matter of isospin asymmetry 6 = {pn — Pp) / {pn + Pp) can be written as 

e(p,5) = e(p,0) + E,,^(p)52 (1) 

within the parabolic approximation (see e.g., [4]), where e(p, 0) is the energy per nucleon 
in isospin symmetric nuclear matter and Egymip) is the symmetry energy at density p. 
To study the Egymip) has been a longstanding goal of extensive research with various 
microscopic and/or phenomenological models over the last few decades. The predicted 
results on the density dependence, especially at high densities, are extremely diverse and 
often contradictory. Theoretical results can be roughly classified into two groups, i.e., a 
group where the Egymip) rises monotonously and one in which it falls with the increasing 
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density above about twice the normal nuclear matter density. Above 3po, even the sign 
of the symmetry energy is unclear. The fundamental cause of the extremely uncertain 
HD behaviour of Esym{p) is the complete lack of terrestrial laboratory data to constrain 
directly the model predictions. Recently, several promising observables were proposed to 
extract the HD behavior of the symmetry energy using high energy radioactive beams [5]. 
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Figure 1. Upper window: the isospin asymmetry-density correlations at t=20 fm/c and 
Ebeam/^ — 400 MeV in the central ^^'^Sn -|-^^^ Sn reaction with the nuclear symmetry 
energy E^y^^ and E'^y^, respectively. Lower window: the same correlation as in the upper 
window but at 10 fm/c and Ebeam/^ = 2 GeV/nucleon. The corresponding correlation in 
neutron stars is shown in the insert. 



2. CORRELATION BETWEEN DENSITY AND ISOSPIN ASYMMETRY 
IN NEUTRON STARS AND HEAVY-ION COLLISIONS 

In a neutron star, several of its main properties arc affected significantly by the cor- 
relation between the density and the proton fraction xp at /3 equilibrium. The latter is 
determined by [6] 

hc{37T^pXfs)'/^ = 4E,y„,(p)(l - 2xfs). (2) 

The equilibrium proton fraction is therefore entirely determined by the Esym{p)- We use 
the following two representatives of the symmetry energy 

3 — -u 

Etymip) = Esym{po)u and Ely^{p) = Esym{po)u ■ , (3) 
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where u = p/po- The values of = 1 — 2x^3 corresponding to these two forms of 
symmetry energy are shown in the insert of Fig. 1. With the E^y^{p), the 5^ is 1 for 
p/ Po ^ 3, indicating that the neutron star has become a pure neutron matter at these high 
densities. On the contrary, with the Egy^{p), the neutron star becomes more proton- rich 
as the density increases. To see the connection between the physics governing properties 
of both neutron stars and heavy-ion colhsions, we show in Fig. 1 the correlation between 
density p and the isospin asymmetry Suke (effective isospin asymmetry taking into account 
different charge states of A(1232) resonances) over the entire reaction volume at the time 
of about the maximum compression in the central ^^'^Sn+^'^'^ Sn reactions with E/A=400 
(upper window) and 2000 (bottom window) MeV/nucleon, respectively. These results 
were obtained using an isospin-dependent hadronic transport model [7]. The overall rise 
of 6 at low densities is mainly due to the neutron skins of the colliding nuclei and the 
distillated neutrons. Effects due to the different symmetry energies are more clearly 
revealed especially at high densities. An astonishing similarity is seen in the resultant 5—p 
correlations for the neutron star and the heavy-ion collision. In both cases, the symmetry 
energy E\y^^{p) makes the HD nuclear matter more neutron-rich than the E'^y„^(p) and 
the effect grows with the increasing density. Of course, this is no surprise since the 
same underlying nuclear EOS is at work in both cases. The decreasing E^y^^p) above 
— l-5po makes it more energetically favorable to have the denser region more neutron- 
rich in both neutron stars and heavy-ion collisions. 

3. n-/n+ PROBE OF THE E,j,^(p) 

It is well known that the 7r~ /tt'^ ratio in heavy- ion collisions depends strongly on the 
isospin asymmetry of the reaction system. It is also qualitatively easy to understand why 
this dependence can be used to extract crucial information about the EOS of neutron- 
rich matter. On one hand, within the A resonance model for pion production from first- 
chance independent nucleon-nucleon collisions [15], the primordial 7r~/7r+ ratio is (5iV^ -|- 
NZ)/{5Z'^ + NZ) [N/Zy. It is thus a direct measure of the isospin asymmetry 
{N/Z)dense of the dcusc matter in the participant region of heavy-ion collisions. It was 
shown that the {N/Z)fiense is uniquely determined by the high density behaviour of the 
nuclear symmetry energy[5]. Therefore, the 7r~/7r"'" ratio can be used to probe sensitively 
the EOS of neutron-rich matter. On the other hand, within the statistical model for 
pion production[16], the tt^/tt''" ratio is proportional to exp [(//„ — pLp)/T], where T is the 
temperature, ^„ and pp are the chemical potentials of neutrons and protons, respectively. 
At modestly high temperatures (T > 4 MeV), the difference in the neutron and proton 
chemical potentials can be written as[17] 

Pn ^ Pp — ^asy ~ ^asy ~ ^Coulomb + T 

where Vcouiomb is the Coulomb potential for protons. At is the thermal wavelength of a 
nucleon and h'^s are the inversion coefficients of the Fermi distribution function. The 
difference in neutron and proton symmetry potentials V^^y — V^^y = 2vasy{p)S, where the 
function Vasy{p) is completely determined by the density-dependence of the symmetry 
energy [5]. It is seen that the kinetic part of the difference pn — Pp relates directly to the 
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isospin asymmetry Pn/ Pp or pn — Pp. Thus within the statistical model too, the tt~ /tt^ 
ratio is sensitive to the {N/Z)(iense- Moreover, the value of tt^/tt^ ratio is affected by the 
competition of the symmetry and Coulomb potentials which all depend on the isospin 
asymmetry of the reaction system. At relatively low temperatures there is a good chance 
to extract crucial information about the symmetry potential. 




t (fm/c) 



Figure 2. The evolution of neutron/proton ratio of dense region (left) and tt /tt^ ratio 
(right) in the central ^'"'^Sn Sn reactions at E^^^^^j A = 200,400 and 1000 MeV. 



Shown in Fig. 2 are the isospin asymmetry of the high density region {n/p)p>pQ (left) 
and the {'k' /'K'^)iike ratio (right) 

as a function of time for the central ^^'^Sn +^^^ Sn reaction at beam energies from 200 to 
1000 MeV/nucleon. This ratio naturally becomes the final vr"/??"'" ratio at the freeze-out 
when the reaction time t is much longer than the lifetime of the delta resonance ta- The 
{tt~ /TT^)iike ratio is rather high in the early stage of the reaction because of the large 
numbers of neutron-neutron scatterings near the surfaces where the neutron skins of the 
colliding nuclei overlap. It is seen that a variation of about 30% in the {n/p)p/>pQ due 
to the different Egymip) results in about 15% change in the final t[~ /t^^ ratio. It has 
thus an appreciable response factor of about 0.5 to the variation of HD n/p ratio and is 
approximately beam energy independent. Therfore, the {t^~ /'n^)iike ratio is a sensitive 
probe of the HD behaviour of nuclear symmetry energy. 




Figure 3. The neutron-proton differential collective flow in the mid-central ^^'^Sn -\- Sn 
reactions at Ebeam/^ = 400 MeV (upper window) and the excitation function of flow 
parameter (lower window) with the nuclear symmetry energy and , respectively. 



4. NEUTRON-PROTON DIFFERENTIAL FLOW PROBE OF THE E,ym{p) 
The neutron-proton differential collective flow is measured by [21] 

I N{y) 

where N{y) is the total number of free nucleons at the rapidity y, p^, is the transverse mo- 
mentum of particle i in the reaction plane, and is and —1 for neutrons and protons, 
respectively. The free nucleons are identified as those having local nucleon densities less 
than l/8po- The Fnp{y) combines constructively the in- plane transverse momenta gen- 
erated by the isovector potentials while reducing significantly influences of the isoscalar 
potentials of both neutrons and protons. Thus, it can reveal more directly the HD be- 
haviour of Esym{p) in high energy heavy- ion collisions. A typical analysis for mid-central 
^^"^Sn -|-^^^ Sn reactions at 400 MeV/nucleon is made in the upper window of Fig. 3. A 
clear signature of the HD behaviour of Egymip) appears at both forward and backward 
rapidities. The slope dFnp/d{ycm/ybeam) at mid-rapidity is different by about a factor 
of 2 in reactions at Ebeam > 200 McV/nucleon as shown in the lower window of Fig. 3. 
This large effect can be easily observed by using available detectors at several heavy-ion 
facihties in the world. 
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5. SUMMARY 

In summary, the HD behaviour of nuclear symmetry energy has been puzzhng physicists 
for decades. In this work, central collisions induced by high energy radioactive beams are 
proposed as a novel means to solve this longstanding problem. Within an isospin de- 
pendent hadronic transport model using two representative density dependent symmetry 
energy functions predicted by many body theories, several experimental probes for the HD 
symmetry energy are studied. Among them, the 7r~/7r+ ratio, neutron-proton differential 
collective flow and its excitation function are found most promising. Future comparisons 
between the experimental data and model calculations will constrain stringently the HD 
behaviour of nuclear symmetry energy. 

This work was supported in part by the National Science Foundation under grant No. 
PHY-0088934. 
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